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The Kraichnan flow provides an example of a random dynamical system accessible to an exact 
analysis. We study the evolution of the infinitesimal separation between two Lagrangian trajectories 
of the flow. Its long-time asymptotics is reflected in the large deviation regime of the statistics of 
stretching exponents. Whereas in the flow that is isotropic at small scales the distribution of 
such multiplicative large deviations is Gaussian, this does not have to be the case in the presence 
of an anisotropy. We analyze in detail the flow in a two-dimensional periodic square where the 
anisotropy generally persists at small scales. The calculation of the large deviation rate function of 
the stretching exponents reduces in this case to the study of the ground state energy of an integrable 
periodic Schrodinger operator of the Lame type. The underlying integrability permits to explicitly 
C3 . exhibit the non-Gaussianity of the multiplicative large deviations and to analyze the time-scales 

^ ' at which the large deviation regime sets in. In particular, we indicate how the divergence of some 

of those time scales when the two Lyapunov exponents become close allows a discontinuity of the 
, large deviation rate function in the parameters of the flow. The analysis of the two-dimensional 

anisotropic flow permits to identify the general scenario for the appearance of multiplicative large 
deviations together with the restrictions on its applicability. 
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I. INTRODUCTION 



The Kraichnan random ensemble of velocities [22] has been extensively used to model various phenomena 
related to turbulent transport both in the inertial interval of scales that develops at high Reynolds numbers 
and at moderate Reynolds numbers where the viscosity effects play an important role [16]. The passive 
transport of scalar or vector quantities in a velocity field is governed by the Lagrangian flow describing the 
evolution of the trajectories of fluid particles. From the mathematical point of view, such flow provides an 
example of a random dynamical system [1] which in the case of the Kraichnan velocities is described by a 
' stochastic differential equation. For the Kraichnan flows corresponding to moderate Reynolds numbers, the 
, methods borrowed from the theory of random dynamical systems or stochastic differential equations appear 

• to provide important information about the transport properties of the flows. To start with, the values 
of the Lyapunov exponents of the flow, whose existence is asserted by the multiplicative ergodic theorem, 
allow to decide whether the flow is chaotic (positive top Lyapunov exponent) or not, leading to different 

C I directions of the cascades of passively advected scalars [12]. More detailed information about the transport 
I properties of the flow may be extracted from the knowledge of the fluctuations of the exponential stretching 
C ' rates around their limiting long-time values equal to the Lyapunov exponents. In the generic case where all 
^ , the Lyapunov exponents are different, the statistics of the stretching exponents may be expected to exhibit 
k>( ' long but finite times a large deviation regime captured quantitatively by a single function of the vector 

rS [ of the stretching rates. Since the existence of such multiplicative large deviation regime is not assured by 

• general mathematical theorems, see [5] for partial results, it is interesting to have at our disposal models 
" " ' where it may be established and studied in detail. 

One such example that has been known for some time is the homogeneous isotropic Kraichnan flow. The 
corresponding stochastic differential equation has been studied in the mathematical literature in the eighties 
and nineties of the last century. In particular, the Lyapunov exponents have been found in [24] and [4] but 
the large deviations have been studied only separately for the top stretching exponent or for the sum of 
them [5]. With the regain of interest of physicists in the Kraichnan model in the mid- nineties, the same 
stochastic equation resurfaced as the model for the Lagrangian flow at moderate Reynolds numbers with the 
motivations, the accents and the language proper to the turbulence theory community [10]. In particular, 
it was realized that many properties of the turbulent transport require more information about the flow 
than the spectrum of Lyapunov exponents and may be expressed in terms of the rate function of the large 
deviations of the stretching exponents. Those include the rate of decay of the moments of advected scalar 
[3] or of growth of those of the magnetic field [11] or, in compressible flows, of the density fluctuations [2], 
the multifractal properties of long-time density concentrations [6] and the threshold for the onset of the drag 
reduction in polymer solutions [9]. In the homogeneous and isotropic Kraichan flow, the large deviation 
regime of the stretching exponents is Gaussian and the corresponding rate function is a quadratic polynomial 
[2, 3]. Its simple form has permitted to extract analytic answers for many characteristics features of passive 
advection in such flows [16]. 

The simplicity of the multiplicative large deviation regime in the homogeneous isotropic Kraichnan flow 
is due to the decoupling of the dynamics of the stretching exponents from that of the cigcn-directions for 
stretching and contraction [3] . As a result also the exact distribution of the stretching exponents may be found 



analytically in this case as it appears to bo related to the heat kernel of the integrable quantum Calogero- 
Sutherland Hamiltonian for particles on the line interacting with the attractive pair potential proportional to 
the function sinh~^ of the inter-particle distance [7, 18]. In the present paper wc analyze the two-dimensional 
Kraichnan flow in a periodic square often used in numerical simulations. The large scale anisotropy due to 
the shape of the flow volume generically persists on small scales inducing isotropy breaking terms in the 
distribution of strain that drives the evolution of the stretching exponents. Due to the presence of such terms, 
the stretching exponents dynamics docs not decouple anymore from that of the (unstable) cigcn-directions. 
The Lyapunov exponents may nevertheless be still computed analytically and their difference expressed in 
terms of elliptic integrals. The distribution of the sum of the stretching exponents is still Gaussian for all times 
(this is a general fact for the homogeneous Kraichnan flows). As for the rate function of the large deviations of 
the difference of the stretching exponents, its calculation may be reduced to that of the ground-state energy 
of the integrable one- dimensional periodic quantum Lame operator. For general values of the coupling 
constant, the eigenvalues of the latter may be found by numerical diagonalization of infinite tridiagonal 
matrices [14]. Those matrices reduce to finite ones at integer values of the coupling and for the lowest 
eigenvalues. Alternatively, the ground state energy of the Lame operator may be found by direct numerical 
integration of the eigenvalue equation. Both approaches permit to obtain the large deviation rate function 
for the stretching exponents that turns out to be non-quadratic although with quadratic asymptotes. The 
analysis of the spectral gap of the Lame operator permits to assess the time scales at which the multiplicative 
large deviation regime sets in. In particular, the divergence of the time scales relative to the multiplicative 
central-limit regime when the difference of the Lyapunov exponents tends to zero accompanies the observed 
discontinuity of the large deviation rate function in the anisotropy parameter at the point where the Lyapunov 
exponents coincide. 

The plan of the paper is as follows. In Sect. II, we discuss the relations between the Lagrangian flow and 

random dynamical systems introducing the concepts of the natural invariant measure and of the tangent pro- 
cess and stating two different definitions of the stretching exponents. In Sect. Ill we introduce the Kraichnan 
ensemble of velocities and discuss how the general concepts considered before simplify in the homogeneous 
Kraichnan velocities. We recall briefly the results about the statistics of the stretching exponents in the 
isotropic version of the Kraichnan model. Sect. IV is the core of the present paper. We discuss there the 
Kraichnan flow in a periodic square, the persistence of anisotropy at small distances, the calculation of the 
Lyapunov exponents and, finally, the large deviations for the stretching exponents, their relation to the Lame 
equation, their non-Gaussianity, and their discontinuity in the anisotropy parameters. The last section col- 
lects our conclusions. Wc believe that the simple model considered here allows to identify a general scenario 
for the occurrence of the multiplicative large deviations when all the Lyapunov exponents are different and 
to understand the mechanism of its failure when some of the Lyapunov exponents get close. 
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II. LAGRANGIAN FLOW AS A RANDOM DYNAMICAL SYSTEM 



A. Random dynamical systems 



We shall start by describing the Lagrangian flow in the language of random dynamical systems. Let us 
consider an ensemble of velocities uf{r) in a bounded region V of the (i-dimcnsional space. Here a; is a 
random parameter belonging to a probability space fl equipped with a probability measure P{diu). Q may 
be taken as the space of velocity realizations. We shall assume the stationarity of the velocity ensemble, 
i.e. the existence of a 1-parameter group of measure preserving transformations a; i— > of f2 such that 
— '^'t+si''')- The Lagrangian flow describing the trajectories of tracer particles carried by the fluid is 
defined by the ordinary differential equation 

^ = <(«). (n.i) 

Under simple regularity assumptions including the spatial smoothness of the velocities it^(r), the solutions 
R^(r) of Eq. (II. 1) parametrized by their time zero position r define a family of random smooth maps 
of the region V such that R^(t") = $"(r) with the composition rule 



(11.2) 



In particular, one obtains a l-parameter group of transformations 
of the product space VxD, which reaUzes the flow dynamics. 



(11.3) 



B. Natural invciriant mesisure 



Note that ^'^^{r) is the time zero position of the solution that at time s passes through r. Suppose that, 
for continuous functions f on V, the limit 

i^J fi'^-si^))dr =: J f{r)^^{dr) (II.4) 

V V 

exists for almost all w and defines a random family of probability measures /x'^ {dr) on V. Note that due 
to the composition rule (11.2), 

j f{^^{r))ii^{dr) = I f{r)f,^^{dr). (II.5) 

V V 

The measures ij!^{dr) describe the distribution of the time zero positions of the Lagrangian trajectories whose 
initial points were uniformly seeded in the far past. For incompressible velocities the uniform distribution 
is conserved by the flow so that ^,"{dr) = j^dr but in the presence of compressibility, the Lagrangian 
trajectories develop preferential concentrations and the measures ^'^{dr) tend to be singular and supported 
by lower dimensional random attractors. 

From the random measures ii'^{dr) one may synthesize a measure M{dr, dco) on the product space VxCl 
defined by the relation 



J f{r,iv)M{dr,du) = j f{r,u) ii^{dr) 



vxsn V 



where ( • • ) denotes the expectation w.r.t. the probability measure P{duj). Note that due to the property 
(11.5), 



j f{^t{r),^t)M{dr,du,) = J f{r,uj)M{dr,du) 



VxQ VxQ 



so that the measure M(dr, dio) is invariant under the 1-paramctcr group (II. 3) of dynamical transformations. 
We shall call M(dr,dio) the natural invariant measure of the random dynamical system (II. 1). Below, 
we shall assume that M{dr,dij) is ergodic with respect to the group (II. 3), i.e. that functions invariant 
under the dynamics are constant M-almost everywhere. 



C. Tangent process 

Muc'li information about the Lagrangian flow may be extracted by looking at the evolution of the separation 
between two inflnitesimally close trajectories. Consider the Jacobi matrix V$^(r) = W^{r) with the entries 

W) = Vj^fir). 

The matrix Wj^(r) propagates the infinitesimal separations: 

5R'{t;r) = W^{r)5r. 

Note that Wo'(r) = Id and that W^{r) satisfies the linear diff'erential equation 

^ = Sfw (n.6) 

with {S^Yj = \7 jU^^{R^{r)) equal to the matrix elements of the strain along the Lagrangian trajectory. 
We shall call W^{r) the tangent process. In studying it below, Eq. (II.6) will play a crucial role. Note 
that the composition rule (II. 2) implies that 

W:+t{r) = W^r«(r))W^-(r). 



In particular, 



W^- (r) = Wr'i^'^-tir))-' 



(11.7) 



We shall be interested in the statistical properties of the d x d matrices W^{r) for fixed but large t 
with (r, w) sampled according to the natural invariant measure M{dr, dco). As for any real invertible dx d 
matrix, one may decompose 

W = O' diag[e^i , . . . , e'''*] O (II.8) 

with a diagonal positive definite matrix sandwiched in between orthogonal ones. One may demand that the 
stretching exponents pi = pff{r) given by half the logarithm of the eigenvalues of the matrices In VF-^VF 
and InW^VF^, be ordered so that pi > ■ ■ ■ > Pd- They carry an important part of the information about 
the tangent process. The joint probability distribution function (PDF) of the time t stretching exponents is 
given by the formula: 

^ d 

Pt{p) = j ^6{pi-pUr))M{dr,dw). 

VxQ. 

Note the relation 

Pt-tir) = ^pT^+iS'^-S.r)) 

between the forward and the backward exponents that follows from Eq. (II. 7). The invariance of the natural 
measure under the 1-parameter dynamics (II. 3) implies then that 

P_t(pi,...,Pd) = Pt(-Pd, .•.,-Pi). 



D. Multiplicative ergodic theorem and multiplicative large deviations 

The main general result about dynamical systems, the multiplicative ergodic theorem of Oseledec [25], see 
also [26], states that under mild assumptions, the limits 

lim - \n{W'^W)'^{r) =: A^(r) 

t — ^iboo t 

exist for M- almost all {rTuo). Besides, 

A^(r) = 0^(r) diag[Ai,..., Ad] 0^(r)^ (II.9) 

where Ai > • • • > A^ are the Lyapunov exponents that, due to the ergodicity assumption, are (r, w)- 
independent. In particular, Aj are the limits when t — > ±oo of the ratios p^^{r)/t for M-almost all (r,a;). 

When all Lyapunov exponents are different, one may expect that for large but finite time the distribution 
of the stretching exponents takes for pi > ■ ■ ■ > pd the large deviation form [3, 18]: 

Pt{p) oc e-*^(^i/*'--w/t) (11.10) 

with a convex rate function H attaining its minimal value equal to zero at the vector A of the Lyapunov 
exponents. In particular, if H is regular around A one would obtain, as a corollary, the multiplicative central 
limit result stating that -^{pt ~ ^t) tends when t — > oo to the vector of normal variables with the inverse 
covariance given by the second derivative matrix H"{X). To our knowledge, no general theorems assure 
existence of the multiplicative large deviation regime, see [5] for some partial results. There are, nevertheless, 
examples of random (and non-random [21]) dynamical systems where the relation (11.10) indeed holds. As 
has been already mentioned in Introduction, the examination of one of such examples is the main topic of 
the present paper. 

For time reversible velocity ensembles, i.e. when the velocities Uf{r) and — M"((r) have the same distri- 
bution, the rate function H possesses the Evans-Cohen-Morriss [15] or Gallavotti-Cohen [17] type symmetry: 

d 

H{-pd/t, -pi/t) = H{pi/t, pd/t) - ^ pi/t , (11.11) 

i=l 

see [2, 16]. The relations of this type were discussed in a closely related setup in [8]. 



E. Stretching along the unstable flags 



We shall call a family F = {Ei)f^^ of z-dimensional subspaccs 

{0} c £;i C • • • C Ed-i CI Ed^Ii!^ 

a flag. An example is provided by the flag composed by the subspaces spanned by the first i 
vectors of the canonical basis of R''. Clearly, the group GL{d) acts on the space of flags. The subgroup 
B C GL{d) composed of the upper triangular matrices preserves the flag Fq and the Grassmannian Or of 
all flags may be identified with the homogeneous space GL{d)/B = 0{d)/D where D = 0{d) fl S is the 
subgroup of the diagonal matrices with entries ±1. We shall denote by dF the normalized 0(d)-invariant 
measure of Or. 

The flow on the volume V induces the flow on the product space V x Or defined by 

Mimicking the constructions from Sect. II B of the natural invariant measure, one may define the measures 
da"{dr,dF) on F x Gr by the relation 

.i!'^^ ^\ I fi'^-sir,F))drdF =: J f {r , F) [dr , dF) 

VxOr VxGr 

if the limit exists for continuous functions / for almost all oj. Clearly, the measures a'^ depend only on the 
past velocities and the formula 

j f{^t{r,F))a-{dr,dF) = J f {r , F) {dr , dF) . (11.12) 

VxOr VxGr 

analogous to Eq. (II. 5) holds. Out of the measures cr'^, one may synthesize a measure T,{dr,dF,dLu) on the 
product space V x Or x Cl by the relation 

J f{r,F,u,)J:{dr,dF,div) = J f{r,F,u,)a^{dr,dF)y 

VxOrxQ VxOr 

For F = OF^ with O G 0{d), consider the Iwasawa decomposition of the invertible matrix W = W^{r) O : 

W = O' c\mg[ci\...,c'''']N 

with O' orthogonal and N upper-triangular with units on the diagonal. The exponents rn = ■qf^{r,F) 
do not depend on the freedom in the choice of O . We shall call them the stretching exponents along 
the flag F. With {r,F,uj) sampled w.r.t. the probability measure S(dr, dF, cL;), they become random 

variables. Their joint time t PDF will be denoted Qtiv)- 

When all the Lyapunov exponents are different then the orthogonal matrices (r) in Eq. (II. 9) are 
determined modulo the right multiplication by matrices from D. In particular, 0'^{r) defines a flag F^{r) 
of subspaces Ef{r) = 0"{r)E^ of (less and less) unstable directions of the flow. In this case, 



I f{r, F) a^{dr, dF) = J f{r, F^{r)) iJi^{dr) , 



VxOr 

i.e. the measure a'^{dr,dF) is concentrated in the direction of the Grassmannian Or on the unstable flags. 
The subspaces Ef{r) may be characterized by the property that for ^ e € E^{r) \ E^_i{r), 

hm -\n\\W-{r)e\\ = Xi 

t— > — CX) t 

with the limit testing the far past asymptotics. The unstable flags F'^{r) depend only on the velocities at 
the negative times and are covariant under the 1-parameter dynamics (II. 3): 

W^{r)F'^{r) = F"'($^(r)). 

We shall call the exponents rif^{r) = r]f^{r,F'"'{r)) the stretching exponents along the unstable flags. 

The exponents r]f^{r) are not equal to the stretching exponents pf^ir) introduced previously. In particu- 
lar, they are not necessarily non-increasing with i although again the ratios 'rifi^{r)/t tend for M-almost all 
(r,w) to the ordered Lyapunov exponents Aj when t ^ ±oo. Although the PDFs Qtiv) and Pt{p) are, in 
general, different (in particular, the latter is non-zero only for ordered arguments), we shall see below that 
the large deviation parts of Pt and Qt are closely related. In fact, the stretching exponents rji are more 
natural objects than the exponents pi and, as we shall see in examples, their evolution is often simpler to 
describe. 



III. LAGRANGIAN FLOW IN THE KRAICHNAN MODEL 



A. Kraichnan ensemble of velocities 

The Kraichnan ensemble of d-dimensional velocities u'^{t, r) is the Gaussian random ensemble characterized 
by vanishing mean and time decorrelated covariance: 

In particular, we shall consider a homogeneous ensemble in the rf-dimensional periodic box V of side L with 
the spatial covariance given by the Fourier series: 

D'^{r,r') = J2 K^-pW^ -{'i--pd)k'k^e'''-^''-'''^d{\k\), (III.l) 

as often used in numerical simulations. The spectral function d will be assumed fast decreasing or of compact 
support so that the resulting covariance and, consequently, almost all velocity realizations uf{r) are smooth 
in space. The simplest would be to take d supported only on the modes with |fc| = i.e. on the lowest 
nontrivial ones. The parameter p in Eq. (III.l) is called the compressibility degree. It is equal to the 

ratio of the covariances ((X^ Vj'u'^')^) ^ (X^(Viu'^-')^) and is contained between zero and one. Vanishing p 
corresponds to an incompressible flow and p = 1 to a gradient one. 

The Lagrangian flow in the Kraichnan ensemble of velocities is defined by the ordinary differential equation 
(II. 1). Because of the white-noise temporal behavior of the Kraichnan velocities, Eq. (II. 1) becomes, however, 
a stochastic differential equation and, in line with more standard notations, will be written in the form 

dR = u^{r)dt. 

In principle, it requires a choice of the stochastic convention, like Ito's or Stratonovich's one, but in the case 
in question both choices lead to the same solutions (this is due to the vanishing of V^jD^-'iO, 0)). As before, 
one obtains from the solutions a family of smooth random maps [23]. 

For the Kraichnan model, the convergence (II. 4) takes place in the norm of the Gaussian process. 
Besides, due to the homogeneity of the velocity ensemble, 

I f{r)M{dr,dcj) = l^j f{r)ii^{dr)^ = ±-Jf{r)dr. (III.2) 

VxQ V V 

Similarly, due to the homogeneity, 

J f{r,F)j:{dr,dF,div)^(^ J f {r , F) {dr , dF)^ = J f{r,F)drx{dF) (III.3) 

VxOrxQ VxOr VxOr 

for some probability measure x{dF) on Or. Note that averaging Eq. (11.12) with respect to the probability 
measure P{du) for functions / independent of r, one infers that 

j {f{Wr{ro)F))x{dF) = J [f{W^(r)F))drx{dF)= jf{F)x{dF), (III.4) 

Or VxOr Or 



i.e. that the measure x{dF) is invariant under the process 



t ■ 



B. Tangent process in Kraichnan velocities 

Further simplifications appear in the statistics of the tangent process W^{r). For positive t, W^{r) 
depends only on the velocities at positive times and n'^{r) on the velocities at negative times. The temporal 
decorrelation of the Kraichnan velocities implies then for any function / of invertible dx d matrices and for 
t > the factorization 



/ 



vxn 



f{W,"ir))M{dr,d^) = l^j f{W^{r))ir{dr) 



V V 

where the last but one equaUty follows from the relation (III. 2) and the last one is again due to the homogeneity 

of the velocity ensemble. Wc infer that it is enough to know the distribution of W}^{rci) for one fixed ro. 
This simplifies considerably the analysis of the statistics of the stretching exponents in the Kraichnan model. 

Let us suppose now that / is a function on GL{d) invariant under the right multiplication of its argument 

by diagonal matrices with entries ±1 so that f{W) for W = Wlf{r)0 depends on O only via the flag 
F — OF^ (but is not, in general, a fimction of W^{r)F only), see Sect. HE. Again due to the decorrelation 
of velocities at positive and negative times and the relation (III. 3), 



j f{W)i:{dr,dF,dij) = j f{Wi'{r)0)a^{dr,dF)^ 

yxGrxn VxGr 

= W\ I {f^^t{r)0))drx{dF) = j {W^r^) O)) x{dF) (III.5) 



VxOr Or 



The last relation permits to simplify the analysis of the statistics of the stretching exponents along the 
unstable flags. 



C. Multiplicative stochastic equation for the tangent process 

For fixed ro, the distribution of the tangent process W^{rQ) may be obtained by solving the multiplicative 
stochastic equation 

dW = S^Wdt, (III.6) 

the stochastic version of Eq. (II. 6), with initial condition Wq = Id. The further crucial simplification, due to 
the time decorrelation and spatial homogeneity of the Kraichnan velocity ensemble, is that in Eq. (III.6) one 
may take {S^Yj = Vju'^^ro) instead of {Sf Yj = VjUf\R^{rn)). In otlicr words, the dependence on the 
trajectory Rf{ro) may be dropped from as long as we are interested in the distribution of W^{ro) for 
fixed ro, provided we consider the differential equation (II. 6) with the Ito convention (here the convention 
does matter, see Appendix A in [16]). 

Summarizing, the strain process in Eq. (III. 6) may be taken as the matrix-valued white noise with 
mean zero and the covariance 



iS't)\ {Sf,yi) = 5{t-t')V^.V,nD'^ (0,0) =: 5{t-t')Cl]. (III.7) 
For the spatial velocity covariance D^^{r, r') given by the Fourier series (III.l), 

= E [{I-pW -{i-pd)k^y]kkhd{\k\) 

= 2a 5% + 13 {SiSi + 5\5i) + 7 S'^S^i (III.8) 
where fc* = ki and is equal to 1 if i = j = k = I and to zero otherwise. The compressibility degree 



2q+ l)J + 7 



trS'^^sA 2a + 2j3 + d-i 



The positivity of the covariance imposes the inequalities 

7>|/3|, 4a+((i+2)/3 + 27>d|/3|. 

Wc shall exclude the trivial case a = /? = 7 = 0. The case of vanishing a corresponds to the isotropic 
situation when the distributions of S*" and of OS'^O^ coincide for any orthogonal matrix O. The a-term 
breaks the 0(d)-invariance of the distribution of to the discrete subgroup of the symmetries of a cube. 
It is the source of a small scale anisotropy that occurs generically in the Kraichnan flow in a periodic box. 
Indeed, vanishing of a requires a fine tuning of the spectral density d(|fe|). For example, when d is non-zero 
only for |fe| = ^ then necessarily a 7^ 0. 



D. Generator of the tangent process 



The generator of the process Wt satisfying stochastic equation Eq. (III. 6), i.e. the operator £ such that for 
any regular function / on the group GL{d), 

is given by the formula 

d 

d 

= a ^(f >)2 + i (/3 + 7)f 2 _ j2 ^ 1 _ + i + + 1 ^ (ni.9) 

i=l 

where, for E"^^ denoting the basic matrices with the matrix elements {E'^^\ = S^'^Sji, 

{S/f){W) = ^U„/(e-^W) = -J2wid^.^ f{W) , = J2^/£/ , 

k i,j 
ij i 

Note that S"^ is the quadratic Casimir of gl{d), J"^ the one of so(d) and T) the generator of the dilations. 
Formula (III. 9) goes back to rcf. [27] where it was discussed for the isotropic incompressible case. For any 
values of the parameters, the generator L commutes with the right regular action 

{RumW) = f{WM) 

of GL{d) on functions on itself. 



E. Stretching exponents in the isotropic CEise 

The isotropic case with a = has been treated in ref. [2, 3], see also [7, 27]. Here C commutes also with the 
left action of 0{d) given by 

{Lof){W) = f{0-^W) . 

In particular, C preserves the space of functions invariant under the left and right action of 0{d), i.e. 
functions f{p) that depend on W only through the stretching exponents, sec Eq. (II. 8). In other words, the 
stretching exponents evolve independently of the angles of the 0{d) matrices in the decomposition (II. 8). 
On functions /(p), the generator C reduces to the operator 

d r\0 i~l d r\ n d r\ 

':. = ^(g^.g-(..-.)|-).?(i:|-r-^|:|-. 

The right hand side is the generator of the diffusion process pt satisfying the stochastic differential equations 
[3] 

dpi = i±L\2coth{pi- pj)dt - ^^^^^dt + Citdt (III.IO) 

where Q is the white noise with the covariance 

{CitQt') = [{l3 + l)5i^ + P]Sit-t'). 

The time t PDF Pt{p) may be expressed in terms of the heat kernel of the Calogero-Sutherland Hamiltonian 
[7, 18]. The Lyapunov exponents are given by [4, 24] 

A, = ai(rf_2i + i)-i^±i^ 



and are all different. The large deviation form of Pt{p) is easy to obtain by the following heuristic considera- 
tions [3]. Since \pi — pj\ for i ^ j grows approximately linearly with time, at long times coth(pj — Pj) « ±1 
and the operator Cp should reduce to the asymptotic form 



-1 ^ l—l l—l 



Similarly, the stochastic equation (III. 10) simplifies at long times to 

dpi = Xidt + Citdt 

The long-time asymptotics of pi is now easy to find leading to the large deviation form (11.10) of the PDF 
Pf(p) with the quadratic rate function [2, 3] 

H{p/t) = ^[iZiT-^if - TOT7(E(t-^0)1 (111.11) 

taking its minimal value at p/t = X and possessing the symmetry (11.11). 

Let us turn now to the stretching exponents q'^{r,F) along flags F = OF^. Since rji are defined by 
the Iwasawa decomposition of the matrix W = W^{r) O, functions of r] may be identified with functions 
/ of invariant under the right action of upper-triangular matrices N with units on the diagonal and 
the left action by orthogonal matrices O (such functions are necessarily invariant under the right action by 
diagonal matrices with entries ±1). The calculation of the average of f{W)) with respect to the measure 
T,{dr , dF, du)) is now simplified by Eq. (III. 5). For fixed ro and O G 0{d), the statistics of W^{ro)0 may 
be found by solving the Ito stochastic equation (III. 6) with the initial condition Wq = O. It follows that 

d 



In the isotropic case with a = 0, the generator £ preserves the space of function invariant under the right 
action by upper-triangular matrices N with units on the diagonal and under the left action of orthogonal 
matrices O and reduces on such functions to the operator 



d 



of the form coinciding for all times with the asymptotic form C^^ of Cp. As the result, 

mtij.)0)) = /e--^(0,r7)/(r,)dr, = ^ 7^!L(,v.....^A) dv 

with H given by Eq. (III.ll). In particular, the 0-dependence drops out and the integral over the Grass- 

mannian Gr on the right hand side of Eq. (III. 5) becomes trivial (in fact, in the isotropic case, x{d-F) = dF 
i.e. it is 0(d)-invariant measure on Gr). We infer that the time t PDF Qtif]) of the stretching exponents 
along the (unstable) flags is Gaussian for all times: 

p-tH{r,^/t,...,n''/t) 

Qt{ri) 



J Q-tH{n^/t,...,r,l/t) 



In particular, its large deviation form coincides with that for the stretching exponents pt except that in the 

latter case, it is restricted to the region where pi > ■ ■ ■ > pd- We shall see below that such coincidence of 
the large deviation statistics for p and r) holds in more general situations. 



IV. KRAICHNAN FLOW IN A PERIODIC SQUARE 

A. Generator of the 2d tangent process 

We shall discuss here the statistics of the solutions of the multiplicative stochastic equation (III. 6) for 2x2 
matrices with the covariance of the white noise strain S'^ given by Eqs. (HI. 7) and (III. 8) with d = 2. The 
overall time scale of the strain is set by 



The distribution of is invariant under the 90° rotations and reflections in the coordinate axis. If O2L is 
the rotation by 45°, i.e. 



Ok = 



1 



1 1 
-1 1 



then S'f'^ = 0^S"0^ has the similar distribution to that of but with parameters 

a' = —a , (3' = a + P , 'y' = a + j 

Note that t' = t, 2j' + a' = 2^ + a and that the compressibility degree is the same for both sets of the 
parameters. We shall call the ratio 

K = \a\T = |q;'|t' 

the anisotropy degree. It measures the difference between the covariances of the the processes Sf and 
relative to (tr5'^'^5'^) and it is contained between zero and one. We shall use the dependence on k 



to measure the influence of the anisotropy on the distribution of the tangent process Wt- For Wi solving 
Eq. (III.6), the distribution of W! = OjLWtO^ will coincide with that of the solution of Eq. (III. 6) for the 
primed values of the parameters. As the result, the distribution of the stretching exponents Pt{r) as well 
as that of rif{r) for the two sets of parameters are identical. Below, we shall then restrict ourselves to the 
case with a > 0. 



The parametrization (II. 8) takes in two dimensions the form: 



W = 











cos ip sin ip 
- sin tp cos ijj 



(we may assume that det W = 1). Hence functions of W may be viewed as functions of two angles and the 
stretching exponents and satisfying the relations 



f{(l>,Pl,p2,1p) 



f{(t> + TT,p2,Pl,4> - 
f{(t),Pl,P2,'4' + 27r) 



that permit to restrict the parameters to the region pi > p2 and < 4>,tp < 2tt. For the generator of the 
tangent process given by Eq. (III. 9) one obtains the following complicated expression: 



£ = 



7Q; 



9pi + dp^ 



2 smcj) coth(pi — P2)d(t> — cos (j){dp^ ~ dp2) — smcj) sinh ^(pi — P2) 
+ |(/3 + 7) + + 2 sinh-2(pi - p2) 92+1 sinh-2(pi - pa) 92 
- 2 cosh(pi - P2) sinh~2(pi - pa) d^d^ + coth(pi - P2) {dp^ - dp^) 



pi 



a, 



P2 



i(2a + 3/3 + 7) 



pi 



a, 



P2 



(IV.l) 



£ is self-adjoint in the l? scalar product with the measure e"''^"''^ sinh |pi — p2\d(f>dpidp2dijj. Note that 
C commutes separately with the translation of </> by tt, with permutation of pi and with the arbitrary 
translations of tp. Since at the end we shall be interested in the distribution of the (ordered) stretching 
exponents, we may right away restrict ourselves to the sector of functions that are periodic in (j) of period tt, 
even under the interchange of p' and independent of ijj. Upon setting 



|(Pi +P2) = r, pi- p2 = p 
C reduces in the action on ^-independent functions to the operator Cj. + 

Cr = j(2a + 3/3 + 7)a2 - i(2a + 3/? + 7)a,. 



with 



■'<t>p 



2a 



sin (j) coth pd^ — cos (p dp 



+ i(/3 + 7) [2 92 + 2 sinh-2 p9| + 2 cothpS^ 



(IV.2) 



(IV.3) 



It follows that at all times the processes rt and (pt, 4>t), starting at ro = = po and ^0 = 0, are independent. 



B. Sum of Lyapunov exponents 



The PDF of rt takes for all times the Gaussian large deviation form since Cr has constant coefficients: 



Pt{r) = 



V'7r(2a + 3/3 + 7)* 



-tHr{r/t) 



(IV.4) 



with the quadratic rate function 



Hr{r/t) 



2a + 3/3 + 7 



In particular, when t ^ 00, the PDF of r concentrates at ^ = —5(20: + 3/3 + 7) which gives the half of 
the sum of the Lyapunov exponents: 



A1+A2 = -(2a + 3/3 + 7). 



(IV.5) 



Note that if one normalizes the Lyapunov exponents by multiplying them by the overall time scale r = 
(a + /? + 7)"^ then one obtains the relation 

(Ai+A2)t = -2p (IV.6) 

stating that the normalized sum of the Lyapunov exponents is directly tied to the compressibility degree. 

The Evans-Cohen-Morriss-Gallavotti-Cohen symmetry (11.11) involves here only the large deviations of 
and reduces to the identity 

Hri-r/t) = Hr{r/t) - 2r/t. 
C. Asymptotic form of the (p, </))-process 



It remains to find the large-time form of the joint PDF of (p and p. Unlike in the isotropic situation, the 
evolution of p does not decouple from the angle (j) and both have to be treated at the same time. Recall 
that we may restrict ourselves to the sector of functions of period tt in (/> and even in p so that we may 
consider only p > (i.e. pi > P2) imposing the appropriate boundary condition at p = (recall that C 
acts at smooth functions on GL{d)). At long times, p/t will still concentrate at a single value equal to the 
difference of the Lyapunov exponents. Anticipating the latter to be strictly positive, we infer that at large t 
the process pt will take predominantly large values oc t. Consequently, we should be able to replace coth p 



by 1 in the expression (IV. 3) for the generator 
asymptotic form: 



and drop the term with sinh p reducing Cfj^p to the 



2a 



sin 6 dch 



cos 6 dr. 



+ 2^dl + (/3 + 7) 



+ 



(IV.7) 



Note on the margin that under a similar transformation applied to operator £ of Eq. (IV.l), all terms 
involving the angle ip drop out implying that tp becomes frozen at long times, in agreement with the 
Oseledec theorem. With an appropriate boundary conditions at p = 0, the operator is self-adjoint with 
respect to the scalar product with the measure e^clpdcf). The restriction to p > will not, however, 
effect the large deviation form of the PDF of p. In what follows, we shall then simplify the things considering 
the operator as acting on functions of p defined on the whole line, restricting p to positive values only 
at the very end. 

Although £J^* still does not preserve the subspacc of functions that depend only p, it does preserve the 
one of functions that depend only on <p. Hence the evolution of becomes at long times independent of 
that of p although the opposite is not true. More specifically, the angle ^ undergoes at long time a diffusion 
on the circle of circumference tt whose generator 



= 2a sm(hd, 



+ 275^ 



(IV.8) 



is obtained by restricting to functions independent of p. Such diffusion process converges exponentially 
fast (we shall find the rate of the exponential convergence below) to a stationary state. The stationary density 



X{<f>) = 



f 




7 + a sin 



7 + a sin (f 



-1/2 



-1/2 



(IV.9) 



is the unique positive normalized solution of period tt of the equation 

2 



^r^x(0) = ( 



2a 



drh sine 



27 5i x('/') = 0. 



D. Difference of Lyapunov exponents 

The difference of the Lyapunov exponents may be found now by applying the following strategy. Suppose 
that there exists a smooth periodic function /(0) of period tt such that 



It follows that 



dt 



A 



dt 



A = const. 



(IV. 10) 



A 



with the exponentially fast convergence. We infer then that A is equal to the mean asymptotic rate of change 
of p that, in turn, is equal to the difference of the Lyapunov exponents. Identity (IV. 10) may be rewritten 
in somewhat more explicit form as 

2a sm^(j> + /? + 7 + f{'j>) = A. 
Integrating the latter equality against xi^)^ '^^ obtain the relation 



A = 



/3 + 7 + 2a sin (p 



Xif) del) 



(IV.ll) 



that fixes the value of A. It is easy to see that Eq. (IV.ll) is also a sufficient condition for the existence of 
the function f{<f)) satisfying condition (IV. 10). Substituting the explicit expression (IV.9) for x{4')^ we infer 
that 



Ai — A2 



13-1 + 2- 



7 + a sin if 



1/2 



dtp 



I 



7 + a sin (fi 



-1/2 



(IV.12) 



d(p 



Recall that 7 > \f3\ and 2a + 2/3 + 7 > \f3\ and that we have assumed that a > 0. It follows that 
A^ — A^ >/3 + 7>0 and at least one of the last inequalities is sharp unless /3 = 7 = 0. Hence Ai > A2 
if /? + 7 > which is consistent with the assumption that, typically, p becomes large for long times. The 
integrals are given by the elliptic functions K{k) and E{k) with the modulus k = 



7 + a sin (fi 



1/2 



dip 



I 



7 + a sin (p 



-1/2 



= (a + 7) 



dip 



Ejk) 
K{k) 



Fig. 1 depicts the normalized Lyapunov exponents AjT for r = (a + /? + 7) ^ as functions of the anisotropy 
degree k = |a|r, 



Air = -1 + 



3 — 2p+ K 



E{ 



2k 



3-2p+K 



K{ 



2k 



A2r 



2p 



3 — 2p+K. > 
3-2p+K E(\/ 3-2p+^) 



K 



(-\/ 3-2p+k) 



for throe values of the compressibility degree p = 0, p = ^ and p = 1. Note that Air decreases and 
A2r increases with k, at constant compressibility degree p, with the sum of the two fixed to — 2p, see 
Eq. (IV. 6). The incompressible system stays always chaotic (i.e. with positive top Lyapunov exponent) and 
this is also true for sufficiently small compressibility degree. For p slightly below ^, however, an increase 
of K may kill chaos. For p > 5 the system is never chaotic. For p = 1, the tendency of anisotropy to bring 
the Lyapunov exponents closer attains its maximum with the two Lyapunov exponents coinciding for the 
extreme anisotropy when S'^ is a diagonal matrix with independent equally distributed entries representing 
independent stretching and contraction along the coordinate axes (or when it is the 45° rotation of such a 
matrix) . 



(a) p = 



(b) P = I 



(c) P=l 



FIG. 1; Normalized Lyapunov exponents Air and A2T as functions of anisotropy degree k — |a|r for three values 
of the compressibility degree p. 



E. Large deviations for exponents p 



The joint PDF of (f>t and pt takes the form of the heat kernel 



-(0,0;(/),p) 



and for large time t should be well approximated by the modified heat kernel 

Pr{<l>,P) = e*^rp(O,O;0,p). 
In the latter, the p-contribution may be diagonalized by the Fourier transform so that we get 



(IV.13) 



where the integration is over a line Re v = —1/2 parallel to the imaginary axis and 



= 2a 



sin — V cos( 



is a second order differential operator acting on periodic functions of period tt. For Reiy = —1/2, £^ is 
self-adjoint with respect to the scalar product with the measure dcj). As we shall see below, the contour 
C of integration in Eq. (IV.13) may be moved to any line parallel to the imaginary axis. Operator Ci, may 
be viewed as a perturbation of the generator C^" of Eq. (IV. 8) with which it coincides for v — Q. By a 
rescaling, a similarity transform and the elliptic change of variables 



«(0) 



a + 7 



7 + a sin ip 



1/2 



dil) , 



Ly is put into the form of a one-dimensional Schrodinger operator: 



2(a-|-7) 



for the function 



/i(0) 



\ In 



and the attractive potential depicted on Fig. 2, 

7 7 



a + 7 7 + a sin 



7 + a sin 



fc cn (m, fc) 



-fc^ -I- fc'^sn-^(u, fc) , 



(IV.14) 



Q + 7 



. The 



where sn(w, fc) and cn(w, fc) are the Jacobi elliptic function corresponding to the modulus fc = 

Schrodinger operator of Eq. (IV.14) acts on periodic functions of u with the period 2K{k) (that corresponds 
to the period tt in 0). Up to a constant, it is equal to the Lame integrable operator in the Jacobian form 
[14, 29] 



Hy = — — + + l)k'^sn'^{u,k) . 
du'' 



FIG. 2; V(<;/>(u)) as a function of for, from top to bottom for w = 0, = ^ = 0.2, 0.6, 0.9 and 0.99999. 



We thus obtain the relation 



dtj) 27ri 
C 

Note that by the Feynman-Kac formula, for v — i^i +W2 with vi 2 real, the absolute value of the integrand 
on the right hand side is bounded by 

g-i/ip + t(/5+7)[i/i(i/i + l)-i'^]+i/i(h(O)-h(0) ^-^a+-,)t{^--^ + [yi(yi+l)-yl\V(4,(u))^ u{(j))) ''"'^'^^ 

For 7 > 0, the last expression tends to zero when 1/2 ±00 (uniformly on bounded intervals of vi) since 
/3 + 7 > and V is attractive. It follows then that the contour of integration C on the right hand side 
of Eq. (IV. 15) may be, as announced, moved to any line parallel to the imaginary axis. Let us consider the 
spectral decomposition 

00 

^-2i»+,)tn. ^ ^e-2("+^)*^-|M'.,„)(vI/.,„| (IV.16) 

ri=0 

with the eigenvalues E^^n = E^,n{k'^) of Tli, ordered in a non-decreasing way for v real. For large t, the 
dominant contribution comes from the ground state ^'1^,0 of ^i/- In particular, for the vanishing coupling 
constant, *o,o = {2K{k))-^/^, E^fi 0, £'0,1 = i^o,2 = TT^K{k)-^ and e*'"" (0; u(0)) converges at long 
times to {2K{k))^^ with the exponential rate equal to 27r^(a + 'y)K{k)~'^ . Note that this rate goes to zero 
when 7 tends to zero since the half-period K{k) diverges in this limit. 

Insertion of the expansion (IV.16) into the expression (IV. 15) permits to extract the large deviation form 
of the PDF of pt from the ground state contribution: 

P^'{^,p)d^ oc ye-*K/*-(2"+/3+7)-(-+l) + 2(a+7)iJ..o(fe^)] ^ ^-tH.ip/t) (jy^^^) 

C 

with the rate function 

Hp{p/t) = max \vp/t - {2a + (3 + + 1) + 2{a + 7) E^^aik^)] (IV.18) 

defined for p/t > 0. We shall denote by Vmax the value of 1/ where the maximum is attained (the maximum 
over real corresponds to a minimum of the real part along the axis parallel to the imaginary one). Note 
that only the ground state energy of operator Ti^, contributes to the rate function Hp. The ground state 
wave function enters the prefactors in the long time asymptotics of the PDF of p. Let us check that when 
t 00 then p/t concentrates at the value equal to the difference of the Lyapunov exponents as given by 
Eq. (IV. 12). To this end, we must find the minimum of Hp{p/t). The stationarity condition implies the 
equations 



= , p/t ^ {2a + f3 + -r){2iy + 1) - 2{a + 7) d,E,,oik^) 



(IV.19) 



and the minimizing value of p/t is given by the relation 

{plt)min = 2a + /3 + 7 - 2(a + 7)a^£'^,o(fc2)U=o. 
The derivative of the ground state energy may be calculated by the first order perturbation theory: 

2K(k) 

5.£;.,ok=o = (*o,o|fc'sn2(.,fc)|*o,o) = fc' + j V{(t>{u)du 



= 1 - 








where the last equality follows from the elliptic identity 2.584.37 in [20]. We obtain this way the relation 

{p/t)min = /? - 7 + 2(a + 7) 

which agrees with Eq. (IV. 12) for Ai — A2 = A. 

A closed expression for the second derivative of the rate function Hp may be obtained by differentiating 
twice the defining relation (IV. 18): 

^p^pI^^ = d{^) = 2(2a + /? + 7 - (a + 7)92i;.,o(A:2)|,=,_J ^^^"^^^ 

since fmax is related to p/t by the second of Eqs. (IV. 19) whose differentiation leads to the last equality. In 
particular, 



dVn 



1 



^"^^^ d{plt) X 2 (2a + /3 + 7 - (a + 7) d^,E,,oik^)\,=o) ^^^'^^^ 
By the second order perturbation expansion, 

00 

5X,o|.=o = 2d,E,,ol=„ - 2^i;o;i |(^'o,o|fc'sn2(.,fc)|^'o,„)|2 

n=l 



n=l 



1 



with the Fourier coefficients fn{k) = J cos{7rnx)sn^ {K{k)x,k)dx. Hp{X) is equal to the inverse variance 



-1 



of the normal random variable obtained by the central limit lim (pi — A)/ Vt 



t — >oo 



Combining Eqs. (IV. 18) and (IV. 4), one obtains the following form of the large deviations rate functions 
for the stretching exponents: 

^^P^/^'P'/*^ - 4(2a + 3/3 + 7) 



+ max 



-^)-(2a/3 + 7)K'^ + l) + 2(a + 7)^.,o(^) • (IV.23) 



for pi > p2- When a = (the isotropic case) then E^^ = and the large deviation rate function (IV.23) 
reduces to the one given by the two-dimensional version of expression (III. 11). 

F. Large deviations for exponents rj 

Similarly as in the isotropic case, the PDF of the stretching exponents r/t along the unstable flag may be 
obtained from Eq. (III. 5) by considering the Iwasawa decomposition 





" 1 






) To' 









w = \ : I \ J u ; • (iv.24) 



of the solutions of the linear stochastic equation (III. 6) with the initial condition given by a (random) rotation 
matrix. In the parametrization (IV. 24) and upon the substitution +^^2) ~ "i^^ 'ni~^2— the generator 
£ of the tangent process takes the form: 

C = i(2a + /3 + 7)[i52 + 2cos'cbd^+ 2sm^<j)dl + 2sin2 ^e-^-^a^ - 2sin(20)9^9^ 
+2 sin(20) e"''a,,a^ - 4 sin^ <j) e'' d^d^ + 2sin2(/)a^ + sin(2(/)) 9^ - 2 sin(2(/i) e^^'c^^ 
i(/3 + 7)[2sin2 092 _ 2sm'^(Pdl + 2 cos^ (j> e'^"^ d"^ + 2 sm{2(j>) d^d^ ~ 2sm{2(j>) e'^^d^d^ 
-Acos^ (jie-'^d^d^ + 2cos^ (j)d^ - sin(20)a0 + 2 sin(20) e-'^sJ + 27 9^ 
+ \(3dl - i (2a + 3/3 + 7) a,.. 

It is self-adjoint in the scalar product with the measure e~'^^^^d(pdrdridx. In the action on functions 
independent on x, C reduces to the sum + with Cr given by Eq. (IV. 2) and 



2a 



sin (f)d^ - cos (/) 9^ + 27 9? + (/3 + 7) 



dr. 



(IV.25) 



Note that the operator C^^^ has the form identical to the asymptotic form of /I^p, see Eq. (IV. 7). In 
principle, it acts now on functions of with period 47r but it preserves the subspace of functions with period 
TT. The evolution of r decouples from that of (j> and 77 leading to the PDF (IV. 4). 



FIG. 3: Vmax as a function of {r]/t)T for k, = 0.8 and p = 0.9. 



As for the joint time t PDF of cf) and 77, it is related to the heat kernel of by the equality 

TT 

Pt(0,ry) = /e*^*''(0o,O;</>,r7)x('/'o)d0o (IV.26) 



obtained from Eq. (III. 5). Indeed, it follows from the relation relation (III. 4) that the probability measure 
x{dF) on the Grassmannian Or of the 1-dimensional subspaces spanned by vectors (cos^, — sin^) has to 
be proportional to x{M d(j)o where the function x(0) is the stationary density (IV. 9) of the angle 0. Since 
x((/>) is periodic with period tt (as a consequence of the invariance of the law of the tangent process with 
respect to the rotations by 90°), so must be Pt( • , 77) and we could restrict the (?!)o-integration in Eq. (IV.26) 
to the interval [0,7r[. The rest of the analysis of the large deviations of rj runs as for the large deviations 
of p. As the result, the large deviation rate function for rj has the same functional form (IV. 23) as that 
for p except of the absence of the restriction r/i > 772- In particular, the rate function iJ.,, of the difference 
77 = r/i — r/2 of the stretching exponents coincides with the function Hp as defined by Eq. (IV. 18) on the 
whole real line rather than on the half line: 

HJr]/t) = max [7^77/< - {2a + (3 + -i)v{v + 1) + 2{a + 7) E^^q{1?')\ . (IV.27) 

Note that Hn{ri/t) + rj/{2t) is then an even function of -q/t. Indeed, the values of Vmax at which the 
maximum on the right hand side of Eq. (IV.27) is attained for rj/t and for — 77/t are related by the reflection 
around v = —1/2. They are smaller than —1 for rj/t < —A, lie in between —1 and for —A < rj/t < A 
and are positive for rj/t > A. Fig. 3 presents the graph of the maximizing Vmax as a function of {ri/t)T 
for K = 0.8 and p = 0.9. We obtain this way the fluctuation relation that compares the rate function for 
opposite values of rj/t: 



H.,{-rj/t) = H.,{Tj/t) + 7j/t 



It resembles the Evans-Cohen-Morriss-Gallavotti-Cohen relation (11.11) but is different from it. It states, 
in particular, that although i]/t may take negative values, the probability of such events is exponentially 
suppressed for large t. 



G. Properties of the rate function Hn 

The large deviations rate function iJ^ is related to the ground state energy Ei,fi{k'^) of the Lame operator 
Til, for k"^ = by the formula (IV. 27). A lot is known about the eigenvalues and eigenfunctions of Tiu. 
The power series expansions for the eigenfunctions in an appropriate parametrization may be obtained by 
solving a recursion relation [29] or by diagonalizing tridiagonal matrices [14]. For v a positive integer and 
the lowest 2p+1 eigenvalues, the corresponding matrices become finite and one obtains as the eigenfunctions 
the "Lame polynomials" . The Lame operator may be diagonalized (for general quasi- momenta) by the Bethe 
Ansatz [13, 19] that in this case goes back to the work of Hermite [29] in the second half of the nineteenth 
century. A simple Maple program computes E^fi{k'^) [28]. We also used a C program to compute E„,o 
by solving directly the eigenvalue equation. Fig. 4 presents the graph of E^flik"^) as a function of v for 
A;2 = 0.2, 0.4, 0.6, 0.8 and 1. 




FIG. 4: E„,oik'^) as a function of f from bottom to top for = 0.2, 0.4, 0.6, 0.8 and 1. 



For large values of f [28] , 

E^e) = iH-i(i + fc') + o{\vi^). 

This leads to the following large behavior of 7?^: 

Note that the left and right asymptotes are parabolas (displaced with each other along the horizontal axis) 
and that 

i/"(±oo) ^ 



2(2a + /3 + 7) 



in agreement with Eq. (IV. 20) since d^E^^ tends to zero at large [I'l. Recall that we have calculated the 
central-limit inverse variance H!^{X) before, see Eq. (IV. 21), so that around the minimum, 

■i 2a + p + J - [a + -/)d^E^^o{k^)\^=o * 

with d^E^ o given by Eq. (IV. 22). The difference between i?^'(±oo) and i?,"(A) attests to the non-Gaussian 
character of the large deviations of r/. Fig. 5 presents the behavior of ^H^{±oo) and of i_ff,"(A) for 
p = 0, 0.7 and 1 as functions of the anisotropy degree k = ar. 

The first quantity diminishes with growing k from the isotropic value ^ for k = to the extremely 
anisotropic one equal to i for k = 1 whereas the second one increases starting from the same initial value. 
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FIG. 5: From bottom to top: iiJ","(±cx)) (the lower curve) and iil/,"(A) for p — 0, 0.7 and 1 as functions of the 
anisotropy degree k. 

It is plausible that for p = 1, ^H'^{X) diverges in the limit n ^ 1. We infer that the undimensionalized 
central-limit covariance ^H^{X) increases with anisotropy and that the probability of large values of r]/t 
decreases slower than if the rate function iJ^ were quadratic with the H'^ equal to its central-limit value. 
Fig. 6 represents the graph Hrj together with the large value asymptotes and the quadratic approximation 
near the minimum for k = 0.8 and p — 0.9. 




FIG. 6: THn{x/T) as a function of x with the large rj asymptotes on the left figure and the quadratic 

approximation around the minimum (the curve tighter for large values) for k = 0.8 and p = 0.9. 
The vertical lines correspond to r]/t = A. 



H. Time scales of the large deviations regime 

It is interesting to look at the time scales at which the large deviation regime for the difference p or 77 of 
the stretching exponents sets in. There were three approximations involved in reducing the exact PDFs to 
their large deviation form. Let us analyze them one by one. 

The first approximation consisted in replacing the PDF Pt{4>,p) by P^^{(j)^p) involving the asymptotic 
form C^^p of the generator C(f,p. The asymptotic form should set in exponentially fast in time with the rate 
given by the difference A of the Lyapunov exponents. This approximation becomes exact when one analyzes 
the difference r] of the stretching exponents along the unstable flags. 

The second approximation consisted of considering only the contributions of the ground state of TLu to 
the kernel of the operator e"^^""'"^-'*^'' in Eq. (IV. 15) with the contour integral along C = {Rev = Vmax\- 
The contributions of the excited states to that kernel decouple exponentially fast with a rate given by the 
spectral gap of 2(a + ^)'H^^^^ for Vmax equal to the value of v that maximizes the right hand side of 
Eq. (IV. 18). We obtain this way a continuum of time scales that depend on p/t or on rj/t. The plot of 
the gap Ejj^i — iJ^^o of T-C„ as a function of v for fc^ = 0, 0.2, 0.4, 0.6, 0.8 and 1 is given in Fig. 7a. For 




— 1 the gap is expected to vanish for —l<i'<0 and its smaU but positive value on the graph is due to 
the slowdown in the numerical algorithm. Indeed, as noticed before, the gap of Tii, for v — or z/=— 1 
is equal to ■K^K{k)^'^ and tends to zero when fc^ tends to 1 (i.e. when 7 approaches 0). For < — 1 
or 1/ > 0, the limiting value of the gap when fc^ tends to zero should, however, be positive and increasing 
with + 1/2|. Fig. 7b presents the gap as a function of for i' — Q, 1, 2, 3 and 4. Since the interval 
Wmax + 1/2| < 1/2 corresponds to p/t < X or to {ri/tl < A, we infer that the decoupling of the contributions 
of the excited states of Hi, to the PDFs of p or 77 on those intervals takes more and more time when 7^0. 

The third approximation in extracting the large deviation form of Pt{p) or QtiVi) consisted in replacing 
the integral in (IV. 17) along the contour C = {Rev — Vmax\ by its saddle point value. This induces the 
correction to Hp[p/t) whose leading term comes from the one- loop contribution 

lln('" + ^ + ^-(" + ^^^-f-°'— ) = ^ln{H"{X)lH"{p/t)) 
2t V 2a + (3 + ^ - {a + j)d^,E,^ol=o ) 2t ^ ^ ^' 

and similarly for Hri{ri/t). As we noticed at the end of the previous subsection, it is plausible that for p = 1 , 
Hp{\) — H'^{\) diverges in the limit k — » 1 causing the divergence of the last correction. 

I. The case with equal Lyapunov exponents 

At the end, let us consider the special case of the potential flow with maximal anisotropy when p = 1 and 
K = 1, i.e. when /3 = 7 = 0. In this case the strain matrix (5*")^ is diagonal with {Sf)\ — ("^4^)^ — ^'^'^ 

{isr)\{sf,)\) = 2aSit-t')({sr)Msr,)\), {{sr)\isr,r,) = o. 

The solution of the multiplicative Ito stochastic equation with initial condition Wq ~ Id takes the form 

t t 

r fiS'^y^ds-at J {S"^ ds - at-: 

Wt = diag eo , e" 

and the stretching exponents pt are given by the formula 

t t t t 

/ j S^(sY^ds- at j S^(s)\ds- at\ / f Su^is^-^ ds - at fS^(s)\ds-at 

pij — max I e" , e" ) , p2f. — min I e° , e" 

This results in the joint time t PDF 

Pt{pi,P2) = -^e-M'^+'^r e-M'^+'^y 0{p,- p,) 
2'Kat 

or, in terms of r = ^{pi + P2) and p = pi — p2 , 

Pt{pi,P2) = -^e-M^+'')\-U^T e{p). 

iTTat 



Similarly the stretching exponents r/t along the flags F — OF^, are given by the Iwasawa decomposition 
(IV. 24) of the matrices W — WtO with the flags F — OF^ distributed with respect to the measure x{dF) 
on the Grassmannian Gr that is invariant under the process VFt, see Eq. (III. 4). Any such measure has to 
be concentrated on the two flags given by the coordinate axis. On obtains then the formula 

/ S^{s)^-^ ds — at f S^{s)%ds — at 

Vit = e° , V2t = e° 

if of'-' is given by the first coordinate axis or the one with interchanged 77^^ if is given by the second 

coordinate axis. In any case, the joint t PDF of the stretching exponents r] takes the form 

Pt{vi,v.) - -i-e--(^+"re--(^+"r 

47rQft 

or, in terms of r = i(?7i + 772) and 7/ = 771 — 772, 

An at 



The PDF of rt agrees with that of Eq. (IV. 4), i.e. rt/t is the normal variable with mean ~a equal to 
the half of the sum of the Lyapunov exponents and with variance a/t. As for pt/i, it is distributed as 
an absolute value of the centered normal variable with variance Aa/t. In particular, the difference of the 
Lyapunov exponents vanishes and the large deviation rate function for pt is quadratic: 

HMt) = (IV.28) 
Similarly, rj/t is a normal variable with mean zero and variance Aa/t and 

The values of the Lyapunov exponents agree with those given by the limiting values of Eqs. (IV. 5) and 
(IV. 12). Recall however, that the non-Gaussianity of the large deviations, as measured by the difference 
i[i/^'(A) — -H"^'(±oo)], was increasing with the growth of the anisotropy degree k, see Fig. 5. For p = 1, in 
particular, iiJ^'(A) was growing with n whereas iiJ^'(±oo) = 2{k+i) decreased to the value \ for n — 1. 
This seems in contradiction with the results (IV.28) and (IV.29) with the quadratic large deviations rate 
functions for p = 1 and n — \ with Hp{p/t) = H'^{r]/t) — j everywhere and not only at infinity. The 
solution of the puzzle lies in the non-uniformity of the large deviation regime when /3, 7 — > and the two 
Lyapunov exponents tend to each other. As we have noticed in the previous subsection, the time scales at 
which the large deviation regime sets in diverge when 7 — > (and, consequently, (3—^0 and A ^ 0). That 
could explain why the limit of Hp{\) = H'^{X) when 7 is not equal to the value of Hp{0) — H^{0) for 
7 = 0. 




FIG. 8: Illustration of conjectured point-wise convergence of tH,-i(x/t) to a;^/8 (diamonds) for a:: > and to /S — x 
(crosses) for x < 0. The dotted line corresponds to /3 = and 7 = a, the solid one to /3 = and 7 = O.Ola. 



Numerical calculations, see Fig. 8, seem to indicate, however, that Hp still converges point-wise to its form 
for 7 = when 7 — > 0. Such point- wise convergence when 7 — s- cannot take place for the large deviations 



rate function iJ^. Indeed, recall that that for 7 > it is Hj^^ri/t) +ri/{2t) that is an even function of rj/t 
whereas for 7 = 0, the rate function H^{r]/t) is even itself. The point-wise convergence of H^j that is a 
function on the whole real line cannot then hold. Instead, for negative rj/t, the rate function H^{r]/t) should 
converge when 7 ^ to Hr,{r]/t) — rj/t. Let us note that the evenness of H^{r]/t) + ??/(2t) for 7 > is a 
consequence of the relations 

Pt(^,r?) = e*^*''(0,0;ry,^) = e^ e*^*" (^i, r?; 0, 0) = e" e*^*''(^;i, 0; 0, -ry) . (IV.30) 

The first one follows from the self-adjointness of the operator C^.q with respect to the measure d(j)dr] 
on the product of the circle by the real line and the second one from the commutation of with the 
translations of r}. Recall that the PDF of r} is given by the integral 

j e*^i">{<l),7])d^ 


into which the initial and the final angles do not enter in a symmetric way so that the equalities (IV.30) 
do not imply that 

7r 7r 

J e*^'!"' {0,0; (f),p)d^ = y e*^*''(O,O;0,-r/)#. 



Nevertheless, for 7 > the last equality holds if the full PDF of 77 are replaced by its large deviation approx- 
imation, the angular asymmetry showing up only in the prefactors related to the ground state eigenfunctions 
of Hu- On the other hand, for 7 = the angular asymmetry does not decouple from the large deviation 
form of the PDF of 77 and conspires to render the latter even. The lack of point-wise convergence of to 
its form for 7 = is a reflection of the singular behavior of the eigenfunctions of Hv when 7^0. 

To summarize, although the distribution of the stretching exponents rj still exhibits large deviation regime 
when 7 = 0, the corresponding rate function is not equal to the limit of the rate functions for 7 > 
signaling that when two Lyapunov exponents coincide, the occurrence of the multiplicative large deviation 
regime becomes problematic. 



V. CONCLUSIONS 



We have examined in detail the tangent process describing the evolution of infinitesimal separations 
between Lagrangian trajectories in the two-dimensional Kraichnan flow in a periodic square. The process 
is driven by the time decorrelated strain whose distribution is, in general, anisotropic, possessing only 
the symmetry with respect of the 90° rotations and axes reflections. Our interest was concentrated on 
the large deviation regime of the stretching exponents p; or t/^ that appear in the matrix decompositions 
W = O' diag[e^S e''^] O or W — O' diag[e''^ , e''^] A'' with orthogonal matrices 0,0' and upper-triangular 
N. The anisotropy couples the dynamics of the stretching exponents to the evolution of the matrices O', 
in contrast to the situation in the isotropic case where the stochastic dynamics of the stretching exponents 
is decoupled from that of the matrices O' and O or O' and N. The stochastic evolution of the matrices 
O' becomes, however, independent, at least at long times, from that of the stretching exponents, attaining 
exponentially fast a stationary state. The latter feeds to the evolution of the stretching exponents in a steady 
fashion permitting them still to attain the large deviation regime. The large deviation rate function for the 
stretching exponents may be expressed in terms of the ground state energy of an operator on the group of 
orthogonal matrices parametrized by the variables conjugate to the stretching exponents. The contribution 
of the excited states to the PDF of the stretching exponents decoiiples exponentially fast with the rate equal 
to the gap of the angular operator. This scenario for the multiplicative large deviations seems quite general 
whenever the Lyapunov exponents are all different, at least in the homogeneous Kraichnan model. When 
some of the Lyapunov exponents become close, some of the time scales for the appearance of the large 
deviation regime as well as the prefactors multiplying the exponential large deviation PDF may diverge. 

What is special about the flow on the periodic square is that the operator on the orthogonal group in ques- 
tion takes form of the integrable periodic Schrodinger operator of the Lame type facilitating the calculations. 
This simplification due to the hidden integrable structure allowed to obtain closed formulae for the Lyapunov 
exponents in terms of elliptic integrals and to analyze the large deviation rate fimction with precision. The 
results of the analysis show that the anisotropy effects lower the top Lyapunov exponent and increase the 
lower one (relative to an overall inverse time scale), with the sum of the two fixed by the compressibility 
degree of the flow. The sum of the two stretching exponents is normally distributed with the covariance again 
fixed by the compressibility degree. The difference of the stretching exponents, however, exhibits in the pres- 
ence of anisotropy non-Gaussian large deviations. Its central limit covariance grows with anisotropy but the 



quadratic large- value asymptotes of the rate function have the top coeificients that decrease with incrciasing 
anisotropy. This non-Gaussian scenario for the multiplicative large deviations applies, however, only when 
the Lyapunov exponents are different. At the extreme anisotropy, when the strain matrix is diagonal with 
independent equally distributed entries, the two Lyapunov exponents coincide and the large deviations for 
the stretching exponents are Gaussian. We have analyzed in more detail this discontinuous restoration of the 
Gaussianity of large deviations. 

For the Kraichnan flow in a periodic rectangle, the multiplicative large deviations may be analyzed similarly. 
The results will be published elsewhere. It remains an open question whether the Kraichnan flow in a three- 
dimensional periodic box possesses a hidden integrable structure that would permit to extend the analysis of 
the present paper to that case. 
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